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ABSTRACT 

As stars evolve along the Asymptotic Giant Branch, strong winds are driven 
from the outer envelope. These winds form a shell, which may ultimately become 
a planetary nebula. Many planetary nebulae are highly asymmetric, hinting at 
the presence of a binary companion. Some post- Asymptotic Giant Branch objects 
are surrounded by torii of crystalline dust, but there is no generally accepted 
mechanism for annealing the amorphous grains in the wind to crystals. In this 
Letter, we show that the shaping of the wind by a binary companion is likely to 
lead to the formation of crystalline dust in the orbital plane of the binary. 

Subject headings: hydrodynamics - planetary nebulae 



1. Introduction 

During the Asymptotic Giant Branch (AGB) phase, strong winds are driven from 
the outer envelope. After ~ 10^ yrs, the AGB envelope is expelled resulting in a proto- 
white dwarf surrounded by a circumstellar nebula. As the hotter core of the star is un- 
veiled, the nebula ionizes and becomes a planetary nebula (PN). Most PNe are highly 
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asy mmetric, displaying com plex morphological structures such as disks and bipolar jets 



[e.g. 



Balick and Frank 



2002 



and references therein). The engine driving the asymmetry 
is thought to begin during the AGB phase or shortly thereafter in the post-AGB phase. 
A binary companion may be responsible for shaping the ne bula. Recent work suggests 



that most, if not all, PNe have incur red a binary interaction (jSorensen and PoUaccd 12004 
De Marco and MoelEoOsI : ISokeil booel ) . 



The formation of equatorial torii, collimated bipolar jets and circumstellar disks may 
be a consequence of additional energy and angular momentum supplied by the binary com- 
panion either through a common envelope phase or by directly shaping the AGB wind 



( iNordhaus and Blackmaru l2006l : iNordhaus et al.l 120071 : iBlackman et al. 



2001} 



A bin ary can 

also process the wide variety of dust species present in AGB winds (IWatersI 12004 ) . Most 



of the grains are amorphous, with crystalli ne silicates also seen at 10-20% abundance from 
stars with particularly high mass loss rates ( Kemper et al. 2001 : Suhl2002). Some post-AG B 
systems possess a torus of crystalline dust at large radius (iGielen and Van Winckell 120071 ). 



Here, we discuss a method of producing the torus of crystalline dust: annealing in shocks 
induced by a binary companion. A wind bl owing past a binary companion will shock, as in 



Bondi-Hoyle-Ljd;tleton accretion (see, e.g. lEdgarll2004l . for a description of Bondi-Hoyle- 



Lyttleton problem). We will show that the shock temperatures can be sufficient to anneal 
the grains with the for med shock lying in the orbit al plane of the binary. Shock heating has 
been proposed before (INakamoto and Miurall2003l ). but without a possible mechanism. 



2. Numerical Study 



Our code is based on the Flash code of iFryxell et al.l (120001 ). an adaptive mesh re- 
finement (AMR) code based around a Piecewise-Parabolic Method (PPM) hydrodynamics 
solver]^. For this work, we use Flash in cartesian co-ordinates. We have added a simple 
NBODY solver, to model the binary. The wind is modelled by resetting all grid cells within 
a distance r^ind of the primary to a density p^^^^, a temperature T^ind; and a radial velocity 
of f„ind with respect to the primary. The gravitational effect of the bodies is subject to 
softening. The softening length for the primary is less than r^^^^; that o f the second ary is 
chosen to be smaller th an the expected Hoyle-Lyttleton ra dius (see, e.g. lEdgarll2004l ). We 
do not include a jet (cf iGarcia-Arredondo and Frankll2004l ). since it is not relevant to the 
present study. 



The source code is available at 



http : / /flash . uchicago . edu/ 
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Close to each boundary there is a damping region, where the gas density and tempera- 
ture are reduced to their ambient values. We do this to ensure that waves cannot reflect from 
the boundaries. We performed two sets of runs. In the first, we used a 1 Mq primary, with 
Pwind = lO-i^gcm-3, T^.^^ = 103 K, v^,^^ = 3.8 x lO^cms"! and r^,^^ = 2 x lO^^ cm. The 
second set of runs contained a 3 M© primary, with T^^-^^ = 2 x 10^ K, v^^^^ = 6.5 x 10® cm s~^, 
and p^ind and r^-^^ unchanged. This choice of parameters ensures that the wind can escape 
from the system. We refine the grid to ensure that r^-^^^ is always covered by 8 grid cells, 
while the secondary's Hoyle-Lyttleton radius is always covered by 4 grid cells. The gas is 
assumed to be adiabatic, with 7 = 4/3. 



Results 



We shall start by discussing a sample calculation, with a 1 Mq primary, a 0.25 M© 
secondary and an orbital semi-major axis of 6AU. This choice makes the Hoyle-Lyttleton 
radius of the secondary approximately 4.5 x 10^^ cm. Since we require at least four grid cells 
across the Hoyle-Lyttleton radius, our grid resolution is just over 10^^ cm in the vicinity of 
the planet's orbit. 

In Figure [H, we show a volumetric density rendering of the system along two axes. 
Th e spiral structure shown in t he view along the z axis is remarkably similar to that noted 
by iMauron and Huggind (120061 ) in AFGL 3068. The tip of the spiral is attached to the 
secondary, forming a bow shock similar to that of Bondi-Hoyle-Lyttleton flow. The shock 
is drawn into a spiral b y the orbita l motion. The structure is also similar to that noted 
by iTheuns and Jorisseru (119931 ) and iMastrodemos and Morris! (119991 ). Figure [2] shows the 
temperature structure of the spiral. Temperatures in excess of 1000 K are generated in a 
narrow region ~ 10^^ cm thick. Although thin in comparison to the computational volume, 
this region is well resolved by our grid. 

We ran a second numerical experiment, identical except for increasing the softening 
length of the secondary to a value substantially larger than the Hoyle-Lyttleton radius. 
This isolated the effect of the secondary from the orbital motion of the binary. The overall 
density structure was little changed, but the temperatures in the spiral dropped dramati- 
cally. Therefore the orbital 'recoil' of the prim ary i s resp onsible for the macroscopic density 
structures, in support of the 'piston ' model of H However, the local deflection of 

the flow by the secondary drives the shock formation, and associated increase in temperature. 



^Compare especially their Figure 4 to our Figure |T] 
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In a third numerical experiment, we increased the binary separation to 10 AU, but kept 
all the other parameters the same as the first. The flow was similar to the first numerical 
experiment, with a slightly reduced peak temperature in the shock.. 

We also performed a fourth numerical experiment, identical to the first except for a 
100 Mj secondary (making a g = 0.095 binary). The spiral density structure was retained, 
but found that the peak temperature region was slightly lower, and spread over a smaller 
volume. 

We performed a similar set of runs for a system with a 3 M0 primary, and found sim- 
ilar behaviour. Since v^^^^ was increased to ensure that the wind could escape, the shock 
temperatures were higher. 



4. Discussion 

O bservations of post- AGB objec ts reveal crystalline dust to be concentrated in the mid- 



plane (iDeroo and Van Winckell 120071 ) . The presence of crystalline dust has been somewhat 
puzzling, since dust grains which form in AGB winds are expected to be amorphous. Other 
authors have suggested that shock annealing might be responsible for crystallisation, but 
without exploring a mechanism for producing the shocks. Meanwhile, many structures ob- 
served around post- AGB stars are attributed to binary evolution. Our numerical experiments 
suggest that the binary and shock models should be combined: crystalline silicates form by 
annealing in the spiral shock formed as the primary's wind passes the secondary. As the 
outward motion of the spiral shock slows, the arms will merge and appear as a torus. 



4.1. Grain Annealing 

There are a number of important scales which must be assessed to determine whether 
such annealing can occur. First are the shock velocities. Examining the output of our Flash 
runs, we find that the velocities are always ~ lO^cms"^. Pre-shock, the velocities can be 
a factor of a couple higher. Post-shock, the velocities can be a factor of a few lower. In 
our calculations, the velocity is mainly radial at all times, so the crossing time of the high 
temperature region (10^^ cm or so in radial extent in the midplane) is on the order of 10^ s. 

Next is the stopping distance of the grains by the gas. Calculating the distance required 
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for a grain to sweep up its own mass of gas, we find 

where we have assumed a standard dust density of pdust = Igcm^^. For expected gas 
densities and dust sizes, this distance is sufficiently small to ensure that gas and dust are 
dynamically coupled; dust grains will not be blown out by radiation pressure independently 
of the gas. The stopping timescale will be similar to the time required to heat the dust 
grains by collisions. This is given by 

10= ( ...J^V ( ^1 (,,^r « (2) 



10 ^^gcm J \1 fim J V 10^ cms ^ 



Hallenbeck et al.l (120001 ) studied the annealing of silicate grains as a function of tempera- 
ture. The process is very sensitive to temperature (it is controlled by the Boltzmann equation, 
with the rate cx exp {—E /kBT}), but if temperatures exceeding 1067K were reached, the 
amorphous grains would anneal to crystals on times tanneai ~ 280 s. At l ower temperatures 



annealing stalled for around 35 h, after an initial burst of crystallisation. iHarker and Desch 



(I2OO2I ) applied this to spiral shocks in proto-planetary discs induced by gravitational insta- 



bilities, concluding that annealing would be possible. 

The final timescale is that of cooling. We have used a simple adiabatic equation of 
state for the gas. In reality, the gas will be able to cool through a forest of molecular line 
transitions. Computing these accurately is complicated, si nce many chemical r e actio ns can 



occur. To make an estimate, we use the cooling curves of ISchirrmacher et al.l (120031 ). who 



considered the wind of a pulsating star. Using figure 4 of their work, we find 

tcooi=f-^^'l-1.4xl0^s (3) 



2mjy(5rad 

where we have estimated the value of the cooling function for a temperature of HOOK and 
a density of 10~^^gcm~^. The cooling rate, Qrad is roughly proportional to the density. 

What are the implications for grain annealing? First, ^anneal ^'^ short compared to 
the other timescales that, so long as temperatures in excess of 1067 K are reached, we can 
assume that annealing is instantaneous. If we could neglect post-shock cooling of the gas, 
then the crossing time for the whole shock is so long that annealing should occur even if the 
shock temperature is only 1067 K. However, cooling is likely to be significant, so successful 
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annealing requires t^^^^^ < t^^^^. Specifically, we require 

iheat ^ Q 9 ( ^\ I Qrad \ ^ 



if^cooi Vl/^my y^lO^ergg ^s" 

(lO-iSgcm-s) (hook) (i06 cms-0 ^ ^ ^"^^ 

Note that Qrad is approximately proportional to density, so this ratio is fairly insensitive to 
the gas density. Smaller grains will heat (and hence anneal) faster. We calculated t^^^^^ for 
1 /im grains, but this is an upper hmit. We expect the actual grains to be smaller, perhaps 
as small as 0.1 /im, making the heating time an order of magnitude less. Such grains would 
anneal easily - particularly when one recalls that each grain trajectory will encounter the 
spiral shock multiple times. 



4.2. Shock Temperature Scaling 

We now estimate how the shock temperature scales with the system parameters. The 
problem may be split into two parts: the wind from the pr imary may be modelled as a 



spherical Bondi wind (see, e.g. chapter 2 of iFrank et al.ll2002l ). A shock is then induced in 



this wind by the presence of the secondary. Post-shock temperatures of 1067 K will anneal 
the dust grains. 

A most important parameter in a Bondi wind is the sonic point. This occurs at a radius 

7-5 X 10^3 T7^ TT^ cm (5) 



In our winds T ~ 1000 K or less, implying that the sonic point will be comfortably outside 
the orbit of the secondary. Since we construct our initial conditions to be escaping and 
supersonic, it is reasona ble to assurn e that our wind solution should always be supersonic 



(the Type 4 solutions of iFrank et al.l ). 



The Bondi solution must be obtained numerically. To proceed, we shall assume that 
the Mach number of the flow remains constant and that the wind velocity remains close to 
the escape velocity. Both assumptions are correct for highly supersonic Bondi flows. The 
launch Mach number of the wind is then 
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and ^ 

T(r) - ( 2 X 1013 e J 

The flow then shocks due to the presence of the secondary. We can use the Rankine-Hugoniot 
equations to find 

^ ^ (27.M^-(7-l))((7-l).M^ + 2) 

Ti (^ + i)2_^2 ^ ^ 

where Ti and T2 are the pre- and post-shock temperatures respectively, and Ai is the preshock 
Mach number. For a strong shock Ai ^ lina7 = 4/3 gas, Equation [8] becomes T2 ~ 
O.lQM^Ti. 

The Mach number of the shock is not A^iaunch, since we must add in the orbital motion 
of the companion in quadrature to obtain the total relative velocity of the gas. This is 
comparatively straightforward, since the velocity of a body in a circular orbit is a factor 

smaller than the escape velocity. This implies that Ai^ ~ Combining 
Equations O [7] and [8], we see that 

^ 1200 (^X( ) " (^i^) " (9) 

\IMqJ V2xl0i3cmy VlOi^cm/ ^' 

This assumes that the wind velocity is always equal to the escape velocity. Our numerical 
experiments had f^jnd slightly larger than the escape velocity, leaving some 'excess' velocity 
which Equation M does not take into account. Consequently, Equation M is rather more 
sensitive to rorb than in our numerical experiments. We emphasise that this calculation is 
only a rough estimate of the temperatures reached, since we saw a wide range of temperatures 
in our numerical experiments. At the tip of the bow shock the temperature was much higher 
(this material would probably be accreted by the secondary anyway). Moving along the spiral 
arms, we found that the temperatures rapidly dropped to values similar to those predicted 
by Equation O 

Our numerical simulations indicate that the shock temperatures depend weakly on the 
secondary mass - two orders of magnitude in secondary mass lead to ~ 20% change in 
temperature. This is not included in Equation M above, but the simulations are however 
consistent with the much stronger predicted dependence on the primary mass. There are 
also two important scales for the secondary in addition to its own radius, namely the Hoyle- 
Lyttelon radius and the Roche radius. In our simulations and in the analytic scaling above we 
implicitly assume that the Hoyle-Lyttleton radius is smaller than the Roche radius. In cases 
where the reverse applies, we would expect an accretion disk to form close to the secondary. 
The shock structure and temperature may depend somewhat on the ratio of these two radii, 
which we have not pursued in the present work. 
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What limits are appropriate for T2? As we have aheady mentioned, we require T2 > 
1067 K for anneahng to occur. However, it must not be too high, or the grains will vapourise 
- 2000 K is a good upper limit. Equation [H] then provides rough constrains on the systems 
which can produce crystalline dust torii. For a 1 Mq star, we would require 0.6 x 10^^ < 
Torb < 1-1 X 10^^ cm (assuming r^j^d = 2 x 10^^ cm). If the primary were a 3M0 star, the 
limits become 5.4 x 10^^ < rorb < 1 x 10^^ cm (assuming the same r^jnd value). 

5. Conclusion 

We have demonstrated that a binary companion to an AGB star can create a torus of 
crystalline dust. The crystalline dust is formed by the annealing of amorphous grains in the 
spiral shock induced by the companion. Such a torus is likely to be slowly expanding, and 
not in Keplerian orbit around the system. 

Although we are certain that our basic mechanism is robust, we have only made rough 
estimates of cooling and heating. We have demonstrated that the relevant timescales should 
permit annealing, but more work is needed. Future calculations should incorporate gas 
cooling. This is not straightforward, since the relevant temperatures and timescales imply 
non-equilibrium chemistry. Ideally, the dust should also be incorporated into the code as a 
separate, coupled component. 

Crystalline dust torii provide strong evidence for binary interactions in AGB winds. 
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Fig. 1. — Volumetric density renderings of the wind emitted by a 1 M© primary, with a 
0.25 Mq secondary in a 6AU orbit. The view along the z axis is shown on the left, that 
along the x axis on the right. The binary orbits in the z — ^ plane 



- 12 - 




Fig. 2. — Temperature structure of the wind emitted by a 1 Mq primary, with a 0.25 Mq 
secondary in a 6 AU orbit 



